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Abstract 

We propose a new variable, Rj, in order to identify exclusive double-diffractive high 
Et dijet production. The variable Rj is calculated using the transverse energy Et and 
pseudorapidity of the jet with the largest Et- For a purely exclusive event the value 
of Rj — > 1, if we were to neglect hadronisaton and the detector resolution effects. To 
illustrate the expected i2j-distribution we also compute exclusive three-jet production; 
and, moreover, include jet smearing effects. By studying the predictions as a function 
of the size of the rapidity interval, Srj, which allows for additional gluon radiation, one 
can probe the QCD radiation effects which are responsible for the Sudakov suppression of 
the exclusive amplitude. In this way we may check, and improve, the formalism used to 
predict the cross sections of exclusive double-diffractive Higgs boson (and/or other New 
Physics) production. 



1 Introduction 

Diffractive processes offer a unique means to discover new physics at the LHC, see for example, 
m 12 El El- An exciting possibility is to search for Higgs bosons in an exclusive reaction, 
that is pp p + H + p, where the plus signs denote large rapidity gaps. This process allows 
detailed measurements of the Higgs boson properties in an exceptionally clean environment and 
provides a unique signature, especiahy for the MSSM Higgs sector, see |S1E|- In particular, the 
Higgs mass and spin-parity determination can be done irrespective of the decay mode, and these 
studies are at the heart of the recent proposal [7j to complement the central detectors at the LHC 
by forward proton taggers placed far away from the interaction point. However, the expected 



1 



event rate is limited; it is strongly suppressed, in particular by a Sudakov form factor necessary 
to guarantee the exclusive final state, see for instance [HI E] • An analogous Sudakov suppression 
enters the predictions for the exclusive production of dijets, 77, etc. The existing diffractive 
Tevatron data (see, for example, the reviews (TUl HH [HI [IHl [HI [IS] and references therein) are 
not in disagreement with the theoretical expectations for these processes, see ITTIITHI IT^ 1^ . 
However a definitive^ confirmation of the mechanism of central diffractive production is still 
desirable. 

Here we examine in more detail the prediction for the important process of central diffrac- 
tive dijet production at the Tevatron. This process is a valuable luminosity monitor for central 
diffractive Higgs production, and for other exclusive processes which may reveal New Physics, 
at the LHC The corresponding cross section was evaluated to be about 10*^ times larger than 
that for the SM Higgs boson. Thus, in principle, the exclusive production of a pair of high Et 
jets (that is pp p+jj +p in the case of the Tevatron) appears to be an ideal 'standard candle' 
for the Higgs. Note, that the CDF measurements have already started to reach values of the 
invariant mass of the Pomeron-Pomeron system in the SM Higgs mass range. This process is 
important on its own right as a gluon factory. As discussed in ||23H |2| the remarkable purity of 
the diffractively produced di-gluon system would provide a unique environment to study the 
properties of high energy gluon jets. Unfortunately, in the present CDF experimental envi- 
ronment, which does not provide tagging of both forward protons, the separation of exclusive 
events is not completely unambiguous. In particular, in addition to the smearing due to the 
jet-searching algorithm and detector effects (see for example, j^) , there are also hadronization 
and QCD radiative effects, which distort the manifestation of the exclusive di-jet signal, see 
for example |2Hl[2ni- Because the reliability of the predictions for the cross sections of central 
exclusive production of heavy mass objects is so important for the prospects of forward physics 
studies at the LHC, it is pivotal to check (whenever possible) all the important ingredients of 
the perturbative QCD approach derived in [HJ [2] . In this paper we focus on how to expose the 
role of the crucial QCD radiative effects which regulate the amount of the Sudakov suppression. 

Recall, that already in QED, it is well known that we can never observe a pure exclusive 
process. For example, the cross section for e~e^ —>■ fi~fi~^ is exactly zero if we exclude the 
photon radiation and additional lepton-pair production which may accompany such events; for 
a review, see [211 • To determine the cross section we must use the celebrated Bloch-Nordsiek [2Z] 
and Kinoshita-Lee-Nauenberg theorems, and calculate the radiative correction accounting 
for the experimental resolution. In experiments with very good resolution the corrections are 
quite large. 

^The observation of exclusive Xc and 77 events (^^ElEl) by the CDF collaboration has been reported at 
the conferences. These results appear to be consistent with the perturbative QCD expectations ^11221, though 
in reality the scale of the Xc production process is too low to justify the use of the perturbative QCD formalism. 
The Tevatron exclusive 77 data are very important. Here we do not face problems with hadronization or with 
the identification of the jets. However the exclusive cross section is rather small. Future precise measurements in 
the diphoton mass interval 10-20 GeV would allow a significant reduction of the uncertainties in the expectations 
for Higgs production, to the order of 30 — 50%. 



2 



An analogous situation occurs when we consider QCD exclusive processes. Here we will 
apply the Bloch-Nordsieck procedure to exclusive diffractive dijet production. That is we will 
allow for additional gluon radiation in some rapidity interval dr], and study how the cross 
section changes as we change the size of 5ri and the energy fraction which is allowed to radiate 
into 5r]. At present, two extreme mechanisms are used to describe central diffractive dijet 
production. First, the formalism for pure exclusive production |H| has been implemented in 
the ExHuMe Monte Carlo [221 • Second, central inelastic dijet production via the inelastic 
interaction of two soft Pomerons, which results in parton-parton scattering at large Et', this 
process is implemented in the POMWIG Monte Carlo [SO]. The dijet distribution is plotted in 
terms of the variable 

% = M,,/Mx . (1) 

In terms of this variable, the first process corresponds to Rjj = 1, since the mass of the dijet 
system, Mjj, is equal to the mass, Mx, of the whole central system. The second process has 
Rjj < 1 since additional radiation (the fragments of the Pomerons) populate the central region, 
that is Mx > Mjj. 



2 A new signature Rj of exclusive dijet events 

Dijet production, with a rapidity gap on either side, has been measured by the CDF collabo- 
ration, both in Run I 31j and in Run II fTH. IT^ IT^ ITU [T3j . at the Tevatron. However there 
may still be some room for doubt whether exclusive dijet production, pp p + jj + p, has 
been actually observed. As mentioned above, there are various effects which strongly smear 
the Rjj distribution, especially in the absence of double proton tagging. The hope was that 
exclusive events would show up as a peak at Rjj = 1. Unfortunately the Rjj distribution is 
strongly smeared out by QCD bremsstrahlung, hadronization, the jet searching algorithm and 
other experimental effects. For example, it was shown, using the ExHume Monte Carlo |32j . 
that only about 10% of exclusive events with Et > 7 GeV have finally Rjj > 0.8, with the 
CDF cuts used in Run I at the Tevatron. 

To weaken the role of this smearing we propose to measure the dijet distribution in terms 
of a new variable 

Rj = 2Et (cosh r]*)/Mx , (2) 

where only the transverse energy E^ and the rapidity t] of the jet with the largest Et are used 
in the numerator. Here rj* = rj — Ym where Ym is the rapidity of the whole central system^. 
Clearly the jet with the largest Et is less affected by hadronization, final parton radiation etc. 
In particular, final state radiation at the lowest order in ^5 will not affect Rj at all, since it does 
not change the kinematics of the highest Et jet used to evaluate (j2I). So despite the emission of 
an extra jet during the final parton shower, we still have Rj = 1. Thus, to see the role of QCD 

^Note that we systematically neglect the effects arising from the transverse momentum of the dijet system, 
which is very small compared to the Et resolution. 
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Figure 1: Central diffractive dijet production; (a) purely exclusive, (b) via soft Pomeron- 
Pomeron interactions, and (c) with a third jet in a given rapidity interval 5rj. The dashed lines 
represent gluons. 



radiation on the Rj distribution, we only account explicitly for additional gluon radiation in 
the initial state. At leading order, it is sufficient to consider the emission of a third gluon jet, 
as shown in Fig. ^ The reason why it is sufficient to consider only one extra jet, is that the 
effect of the other jets, which, at LO, carry lower energy due to the strong ordering, is almost 
negligible in terms of the Rj distribution. The rapidity Ym is sketched in Fig. |21 In Section 5 
we will compute the exclusive three-jet cross section for different choices of the rapidity interval 
St] containing the jets. 



3 Resume of the calculation of exclusive dijet production 

To compute the Rj distribution we first calculate the cross section of the exclusive dijet pro- 
duction of Fig. HJ We have cTexci = '^o' where jH] 
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The first factor, S"^, is the probablity that the rapidity gaps survive against population by 
secondary hadrons from the underlying event, that is hadrons originating from so/it rescattering. 
It is calculated using a model which embodies all the main features of soft diffraction [33]. It is 
found to be S*^ = 0.026 for pp ^ p + H + p aX the LHC The remaining factor, however, 
may be calculated using perturbative QCD techniques, since the dominant contribution to the 
integral comes from the region Aqqj) -C -C M^. The probability amplitudes, fg, to find the 
appropriate pairs of t-channel gluons {Q, qi) and (Q, ^2), are given by the skewed unintegrated 
gluon densities at a hard scale fi ~ Mh/2. 

Since the momentum fraction x' transfered through the screening gluon Q is much smaller 
than that [x) transfered through the active gluons (x' ~ Qtl \fs <^ a; ~ Miil\fs <^ 1), it 
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Figure 2: The rapidity Yu of the central system. It does not necessarily occur at y = 0. The 
rapidity interval containing the jets is denoted by 5?7, outside of which there is no hadronic 
activity. 



is possible to express fg{x,x' ,Q^, fi"^) in terms of the conventional integrated density g{x). A 
simplified form of this relation is jH] 



fg{X,x',Qf,fi'^) = Rg 



d 



dlnQl 



(4) 



which holds to 10-20% accuracy. The factor Rg accounts for the single logQ^ skewed effect. It 
is found to be about 1.4 at the Tevatron energy and about 1.2 at the energy of the LHC 

Note that the fg's embody a Sudakov suppression factor T, which ensures that the gluon 
does not radiate in the evolution from Qt up to the hard scale /i ~ Mh/2, and so preserves the 
rapidity gaps. The Sudakov factor is jSH E] 



Tg{Qt,fJ') = exp 



271 kf 



l-A 



zPgg{z)dz + / 5Z ^lai^)'^^ 



(5) 



with A = + kt). The square root arises in (0} because the (survival) probability not to 

emit any additional gluons is only relevant to the hard (active) gluon. It is the presence of 
this Sudakov factor which makes the integration in Q infrared stable, and perturbative QCD 
applicable. 

It should be emphasised that the presence of the double logarithmic T-factors is a purely 
classical effect, which was first discussed in 1956 by Sudakov in QED [SH]- There is strong 
bremsstrahlung when two colour charged gluons 'annihilate' into a heavy neutral object and 
the probability not to observe such a bremsstrahlung is given by the Sudakov form factor. 
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Therefore, any model (with perturbative or non-perturbative gluons) must account for the 
Sudakov suppression when producing exclusively a heavy neutral boson via the fusion of two 
coloured/charged particles. 

In fact, the T-factors can be calculated to single log accuracy The coUinear single 
logarithms may be summed up using the DGLAP equation. To account for the 'soft' logarithms 
(corresponding to the emission of low energy gluons) the one-loop virtual correction to the 
gg ^ H vertex was calculated explicitly, and then the scale fi = 0.62 Mh was chosen in 
such a way that eq.® reproduces the result of this explicit calculation It is sufficient to 
calculate just the one-loop correction since it is known that the effect of 'soft' gluon emission 
exponentiates. Thus (0) gives the T- factor to single log accuracy^. 

4 Calculation of exclusive 3-jet production 

Here we consider the emission of a third jet described by the variables x and pt- The variable 
X is the fraction of the momentum of the incoming gluon (denoted by xi in Fig. ^c)) carried 
by the third, relatively soft, jet; that is a; = 1 — x\lx\. The explicit formula for the LO third 
jet radiation can be obtained using the helicity formalism reviewed in Ref . j3H] . We outline the 
calculation in the Appendix, where the general formulae for the exclusive three-jet production 
amplitude are presented; that is, not restricted to LO. In the double logarithm limit, with 
•pt ^ Et and x <^ 1, the exclusive 3-jet cross section is simply the exclusive dijet cross section, 
a'^^^ multiplied by the classical probability for soft gluon emission 



Note the extra factor 1/4, which reflects the suppression of soft gluon emission in comparison 
with the usual classical result given by the expression in brackets. Naively we might expect a 
colour factor iVc, but instead we have This is due to the absence of the colour correlation 

between the left (amplitude M) and the right (amplitude M*) parts of the diagram for the 
cross section, in our case with a colour singlet s-channel state. 

■^Of course, in the case of QCD, the exponentiation of soft emission requires some clarification. Because of 
the non-Abehan structure of QCD, there are indeed some particular cases when the soft-emission factorization 
and Poisson distribution theorems do not hold. This was exemplified, in particular, in Ref. However we are 
interested in a phenomenon of a completely different (classical) nature. In we discussed the NLO correction 
to the double log term caused by the classical current, where the soft gluon radiation exponentiates. This 
accounts for the effect of the energy- and angular-ordered additional soft gluon radiation, which, due to QCD 
coherence, is just part of the cascade generated by the 'primary' gluon. Summation of such soft 'single' logs is 
performed analogously to the DGLAP approach, which results in their exponentiation. This situation is of the 
same nature as the well known Modified Leading Logarithmic Approximation, which, for example, is discussed 
in detail in the book by Dokshitzer et al. |37| . 




(6) 
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If we just keep the coUinear logs with respect to the beam direction, that is we keep the 
condition Qt < Vt ^ Ex, but do not impose x -C 1, then the 3-jet cross section becomes 



97T a^^{E^)\ 1 



dt \ 4:71 Pt X 

where 

a 



AE^ j 2 



3 l^x\\-2EllMl) 



The first term, in the round brackets in (jH)), is the known cross section for the exclusive colour- 
singlet ^f^f-dijet production. The variable t in (|7j) denotes the square of the four momentum 
transferred in this exclusive colour-singlet gg high i?T-dijet process. In other words t is 
measured between the highest Et jet and the incoming gluon which produces the high Et 
dijet system. The last term in round brackets in ((7j) is just the double-log expression for 
the emission of the third jet, see Finally, the factor in square brackets in (jHJ accounts 
for the polarization structure of the 3-jet system. Recall that the exclusive double-diffractive 

kinematics selects events with the same helicities of the incoming gluons, either (++) or ( ), 

that is Jz = 0. The first term, (1 — x)^, corresponds to the helicity of the soft (third) jet being 
equal to the helicities of the incoming gluons, whereas the remaining expression corresponds 
to the third jet having opposite helicity to that of the incoming gluons. In this expression, the 
term proportional to originates from the high E^ dijets having different helicities, whereas 
the factor 1 in the numerator corresponds to the production of two high Et jets with the 
helicities equal to each other. The 1/(1 — x) in the second term reflects the usual (BFKL-like) 
1/z singularity in the Altarelli-Parisi splitting function P{z). 

It is informative to note that the behaviour of all three terms in the square brackets of Eq. 
(JS)), in the x ^ or x — 1 limits, is not accidental. Its physical origin can be understood 
by recalling the celebrated Low soft-bremsstrahlung theorem j^H] (see also [101 141p. Recall, 
that according to the MHV rules (see the Appendix), the only non- vanishing Born 2^2 
amplitudes, Mb, are those which have two positive and two negative helicities. On the other 
hand, the = selection rule requires that the two incoming gluons have the same helicities, 

either (-I--I-) or ( ). According to the Low theorem j^Hj, for radiation of a soft gluon with 

energy fraction z <^ 1, the radiative matrix element Mrad may be expanded in powers of z 

^ oo 

Mrad~-5ZC„^", (9) 



where the first two terms, with coefficients Cq and Ci (which correspond to long-distance 
radiation), can be written in terms of the non-radiative matrix element Mb- 

The application of these classical results is especially transparent when the cross sections are 
integrated over the azimuthal angles. Then the non-radiative process depends only on simple 
variables, such as the centre-of-mass energy ^. In particular, if Mb = 0, the expansion starts 

"^Note that in our case, in the colhnear log approximation, when Qt <^ pt <^ Et, the azimuthal angular 
dependence is practically absent. 
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from the non-universal C2-2^ term, which corresponds to non-classical (short-distance) effects, 
not related to Mb, see lIUl liT]. 

Let us start with the third term in the square brackets in Eq. (jSJ. In this case soft 
radiation should be considered with z = x <^ 1. The corresponding non-radiative matrix 

element vanishes, since its helicity structure is either (+ + H — ) or ( h). Therefore, the 

matrix element squared, jMradP, is proportional to x"^. Keeping in mind the factor x"^ dx/x, 
which arises from phase space, we see that this term is indeed proportional to a;^ dx/x, as it 
appears in Eq. (jHJ. The soft-radiation limit of the first term corresponds to 2; = (1 — x) <^ 1. 
Then the third jet carries the largest momentum, and one of the final jets is very soft. Again, 
the corresponding Born amplitude vanishes due to the MHV rule, and we arrive at the result 
l-^radP ~ (1 — s)^ d{l — x)/{l~x). Finally, the second term, with the factor 1 in the numerator, 
corresponds to the only non- vanishing non-radiative amplitude, either (+ H ) or ( 1-+). 

In the case of the coUinear LO process (i.e. pt <^ Ex), the value of Rj can be calculated as 

^ Vcosh(r7* ± iln(l -x))y ^ ^ 

Here a/1 — x = Mjj/Mx accounts for a lower mass, Mjj, of dijet system in comparison with 
the mass Mx of 3-jet system, whereas the factor in brackets accounts for the corresponding 
shift (by 0.51n(l — x)) of the rapidity of dijet system. The minus sign must be used in (fTUI) 
when the highest Et jet goes in the same (beam or target) hemisphere as the soft (third) jet. 

5 How the third jet affects the distribution in Rj 

With knowledge of the luminosity, Q, and the cross section of the hard subprocess, (jHl), we 
can calculate the cross section of exclusive 3-jet production, and study how this contribution 
looks in terms of the Rj variable. Note that, after the emission of the third jet, the production 
of other soft jets with x' < x practically does not alter the value of Rj. 

In the naive, a la QED, case, this multijet emission cancels a large part of the Sudakov 
T-factor suppression. In other words, it gives an exponent analogous to that in (0), but with 
a positive power. In QCD the situation is more complicated. In the expression for the cross 
section, MM*, the two active t-channel gluons (one in M, the other in M*) are not correlated 
with each other, but form colour singlets, each with the corresponding screening gluon in its 
own amplitude, M or M*, see Fig. |3] The colour decomposition of the t-channel pair of active 
gluons, gg', is given by 

i 

where Ai denotes the colour multiplet of the t-channel gg' system: that is, Ai is the colour 
singlet. As and Ag {Aiq and Aiq) are the asymmetric and symmetric colour octets (decuplets) 
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Figure 3: The cross section, MM*, for exclusive three-jet production, where the active gluons 
are denoted by g and g', see the t-channel decomposition of eq. (fTTj) . The two outside vertical 
lines are the screening gluons; indeed all the lines in the plot denote gluons. The dashed line 
is the third (soft) jet, with kinematic variables x and pt- The colour labels a, b, c, d, e are those 
used in eq. (|T!?|l . 



components, etc. The coefficients c- give the probability to have one or another colour state. 
Thus the probability that the pair of active gluons, gg', forms the corresponding colour multiplet 
is 

Ci = tc'i, (12) 

that is c'i times the statistical weight given by the number i of members of the multiplet. 

If we use the decomposition of the product of two 3-gluon vertices fahefcde over the colour 
projection operators Pj, that is 

I'fabefcde = f 3Pi + ^Pg + \p-^ - P27) , (13) 

V ^ ^ / ab,cd 

then we see that for each t-channel colour multiplet, the probability of soft gluon emission is 
driven by its own colour factor Aj. Namely, we have Ai = A^c = 3 for the singlet. As = 3/2 
for the octets, Aio = for the decuplets and A27 = — 1 for the 27-multiplet. The colour labels 
a, 6, c, d, e are shown in Fig. IHl 

So to compute the Rj distribution we must include the factors arising from including the 
third jet with the corresponding colour charge Aj for each term in the decomposition (jllll . The 
power of the exponent for this real emission has the form of the T for the virtual corrections 
(0) multiplied by the corresponding colour factor Xi/N^.. For instance, for the case when the 
(yf^f'-pair form a singlet, that is for i = 1, we have Xi/N^ = 1. Taking each exponent with its 
weight Ci, we obtain 

r.-. . ^ , exp (I /; P,,,.).. .(.,/2 - W)) , (14) 
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where the scale n = 0.62y M|^/(l — a;) is taken to be the same as in (0) and where the 
coefficients the weightings in the decomposition shown in eq. (jllj) . Unhke eq. 

the z integral is limited by the momentum fraction x carried by the soft third jet; for the case 
of X > 1/2 the upper limit x in the z integral (14) is replaced by 1 — x - two jets cannot carry 
the fraction of an initial momentum greater than 1 (i.e. x + z < 1). Next, we have added the 
^-function, which enables us to vary the size of the Srj interval containing the jets, so that we 
can study the radiation effect in more detail. As a rule, the jet reconstruction is performed in 
some limited rapidity interval, so it is natural to select events where all the jets are emitted 
within the interval Sr] centred at the position of the Mx system (that is in the interval ±Srj/2 
in the frame where Ym = 0, see Fig. while any hadron activity outside the interval 61] is 
forbidden. 

Note that, due to a more complicated colour structure in QCD, even in the double log limit, 
there is no exact cancellation between the real emission ()14j) and the Sudakov T-factor 

To calculate the exclusive cross section for 3-jet production accompanied by the emission 
of softer jets in the rapidity interval Srj, we multiply the exclusive luminosity (jS)) by the cross 
section of the hard (LO 3-jet production) subprocess, (jH)), and by the factor T^"^^^^)^ (El)) to 
account for the allowed radiation of softer gluons. The results are presented in Fig. ^ and Fig. El 
in terms of distributions over the new variable Rj. In order to do this, relation (|T0|) was used 
to transform the distributions over the momentum fraction x carried by the soft gluon, into the 
i?j-distributions presented in the figures. 

To be explicit the procedure is as follows. The Rj distribution is computed using 

^./..|.,.rfr(^),T<-v(§)- (la) 

where the luminosity C is given in Q and the Sudakov factor T^'^^^^^ is given by (fT^ : and where 
r]i and 772 are the rapidities of the high Et jets. We integrate over the kinematic intervals 

Et > -Emin, |?7i,2| < 2.5, Pniin < Pt < Pmax- (16) 

The lower limit of the logarithmic pt integral is given either by the transverse momentum Qt 
in the gluon loop^ or by the allowed rapidity interval 5r], that is Pmin = max {Q^, xM^e"''''/^}. 
The upper limit is of a pure kinematical nature: Pmax = min xMx/2}. If Pmax < Pmin, 
then there is no LO contribution. 

^The simplest example of this lack of cancellation is exclusive Higgs boson production, where already at the 
first as order there is Sudakov suppression 101, while it is impossible to emit only one gluon accompanying the 
Higgs boson from the colourless two gluon state. 

^For Pt < Qt, the destructive interference between emissions from the active gluon xi and from the screening 
gluon (that is, the left gluon in Fig.^a,c)) kills the logarithmic pt integration. Strictly speaking the values of 
Qt in the amplitudes M and M* may be different, but this effect is beyond the LO accuracy of our calculation. 
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Figure 4: The Rj distribution of exclusive two- and three-jet production at the Tevatron. 
Without smearing, exclusive two-jet production would be just a 5-function at Rj = 1. The 
distribution for three-jet production is shown for different choices of the rapidity interval, 5rj, 
containing the jets; these distributions are shown with and without smearing. The highest Et 
jet must have > 20 GeV. 
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Figure 5: The Rj distribution of exclusive two- and three-jet production at the LHC. Without 
smearing, exclusive two-jet production would be just a 5- function at Rj = 1. The distribution 
for three-jet production is shown for different choices of the rapidity interval, Srj, containing 
the jets; these distributions are shown with and without smearing. The highest Et jet must 
have Et > 50 GeV. 
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Next, we have to include the emission of the third (soft) jet in the direction of one or the 
other incoming gluons, that is beam protons. In other words we must sum up the contributions 
with either the plus or minus signs plus in (jl(J|) . Thus, finally, we obtain 

^ = ^ / dEldr^^dr,, E E ^ (^) 
f dQt a;'^^ g2^ /i^)/s(a;2, 4, Ql, /i^) exp(ni)^ln(p^^^/p^i^ 



^ 1 (^^) 



where rii denotes the power in the exponent in t'^'^'^^^) of ()14|1 . The quantity a arising from the 
hard gg ggg subprocess is given by (jH)). Note that the factor dpf/pf in ((Zj) gives rise to the 
logarithm in C in (jl7p . while the factor dx/x goes into {xdRj/dx)~^. Indeed, the value of x 
and the derivative 



dRj Rj 



dx 2(1 — x) 



±tanh ( T]* ± ^ln(l -x) ) - 1 



(18) 



are calculated according (fTUj). Note that since the lower limit, Pmin, of the integration over the 
Pt of the soft jet may depend on the transverse momentum, Qt, in the internal gluon loop, the 
factors exp(ni) and ln(p^ax/Pmin) occur inside the 'luminosity integral'. 

In the computation we have used the partons of Ref . 02] • We neglect hadronization effects, 
and present the parton level results by dashed curves. In terms of the Rj distribution, the 
exclusive dijet contribution occurs as a 5-function, 6{Rj — l), and cannot be shown in the figures. 
However, in any realistic experiment, the distribution is smeared, at least by fluctuations in 
the calorimeter ^. To see the effect of more or less realistic smearing, we assume a Gaussian 
distribution with a typical resolution^ a = 0.6 /y/ Ex in GeV. 

The results obtained, after this smearing of the parton level distributions, are shown by the 
continuous curves in Fig. IH and Fig. El We see that for the case of 6r] < 5 the exclusive dijet 
production still dominates for Rj > 0.7 — 0.8. The perturbative QCD radiation is suppressed 
by the extra coupling a^. However this suppression is partly compensated by the coUinear logs 
and by a large longitudinal phase space, that is by the rapidity interval 6r] allowed for the 
emission of the extra soft jets. Indeed, we see that the cross section grows with 6t], and by 
St] > 10 is close to the saturation curve (denoted "all St]"), which covers the whole interval of 
leading log QCD radiation. 

''If we assume that the two forward protons are tagged, (as is possible, in principle, in DO experiment at the 
Tevatron |43[l44| or at the LHC if the CMS and/or ATLAS detectors are supplemented by the Roman Pots) then 
the mass of the whole system, Mx, can be measured with much better accuracy by the missing mass method. 

'^We thank M.G. Albrow, D. Alton, M. Arneodo, A. Brandt, C. Buttar, R. Harris, C. Royon and K. Terashi 
for discussions on this choice. The resolution a = 0.6/ \/Et in GeV is close to that obtained for the CDF 
detector, namely a = 0.64:/ ^/Et in GeV + 0.028. The resolution of the DO hadron calorimeter is not quite 
so good: a ^ 20% for Et = 20 GeV. Moreover the expected resolution of the CMS hadron calorimeter is 
about twice worse, while the anticipated resolution of the ATLAS detector may be even a bit better: a ^ 
0.5/y/ET in GeV + 0.015. 
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Figure 6: A schematic diagram of the gg 
third jet. 



ggg process. The gluon labelled by e is the (soft) 



Note that in the region Rj < 0.6 — 0.7 the dominant contribution comes from three jet 
emission. Moreover here the results are more weakly dependent on possible smearing. Of 
course, in the region of small Rj there may be other contributions coming from the three- 
or four-jet Mercedes-like configurations^. However these contributions are not expected to be 
large, since in this case as is not compensated by large logs. Another possible contribution 
comes from configurations which look like inelastic dijet production in the collisions of two soft 
Pomerons. Such configurations, corresponding to Fig. ^b), may populate the low Rj region, 
and are beyond the scope of the present analysis. 



6 General use of Rj 

In spite of the fact that the Rj variable was introduced to select exclusive dijets in double- 
diffractive hadron-hadron interactions in which both of the outgoing protons are tagged, a 
similar idea can be used to improve the measurements of the light-cone momentum fraction 
carried by the dijet system in other situations. In particular, to measure the fraction of the 
photon momentum, x^, carried by the high Et dijets in DIS. Note that the final state radiation 
(and hadronisation) affect mainly the energy, and much less the rapidity of the jet. Therefore 
to calculate (or x^- and x~- in the more general case) one can use the Et of the largest Et 
jet together with the rapidity of each jet. 



Appendix: Helicity amplitudes for gg ggg 

Here we outline the formalism used to calculate the gg ggg process shown in Fig. IHl We 
denote the colour indices of the incoming gluons by a,b, and of the outgoing high Et gluons 

^In the Appendix we give the formulae needed to compute exclusive three-jet production in the whole 
kinematical interval, and not just in the domain of the leading coUinear log approximation. 
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by c, d. Finally the colour index of the soft jet is denoted by e. The gg ggg matrix element, 
which depends on the helicities, hi, and the 4-momenta, Pi, of gluons, is given by the so-called 
dual expansion (see [HH] and references therein) 

^hMMe ^p^^ p^^ p^^ p^^ ^ J2 Ti{XaXbXAdXe) m{a, b, c, d, e) , (19) 

where the sum is over the non-cyclic permutations of a, b, c, d, e. The first factor looks as if 
all the gluons were emitted from the quark loop; where Aj are the standard matrices of the 
fundamental representation of SU(3), which are normalised as follows 

Tr(A"A^) = (20) 

[A^A^]=^^,eA^ (21) 

The colour-ordered subamplitudes, m(a, b, c, d, e), are only functions of the kinematical variables 
of the process, i.e. the momenta and the helicities of the gluons. They may be written in terms 
of the products of the Dirac bispinors, that is in terms of the angular (and square) brackets 

{ab) = {p-\pt) = v^2^e^^-, (22) 



Wb] = (pM = V|2paPfc|e^^-, (23) 

where 2[paPb) = Sab is the square of the energy of the corresponding pair. If both 4-momenta 
have positive energy, the phase (pab is given by 

with pf = p'i + Pi, while the phase 0afe can be calculated using the identity Sab = {ab)[ab]. 
Actually the phase (pab is irrelevant in our collinear LO calculations, except for the fact that 
{ab) = —{ba) and [ab] = — [ba]. However to calculate the gg — > ggg amplitude beyond LO, and 
to compute a more precise cross section, based on eqs. P9I25|) . we would have to account for 
the phases. 

Finally, the only non-zero subamplitudes 

m(a, &, c, d, e) = ^9'2'''Tl^rjTS%T\r-\ ^^5) 

{ab) {be) {cd) {de) {ea) 

are those which have two helicities of one sign, with the other three of the opposite sign, 
the so-called Maximal Helicity Violating (MHV) amplitudes. Here g is the QCD coupling 
{as = g'^/47c). In particular, when ha = hb = —1 while h^. = h^ = h^. = +1 the numerator 
(JJ)^ = (a6)^; i.e. / and J are the only two gluons with the same helicities. If we change the 
sign of helicities, then we have simultaneously to replace the {ij) brackets by the [ij] brackets. 
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Note that the colhnear logarithm in the direction of gluon a comes from the factor (ae) (or 
(ea)) in the denominator of (j^ . Thus to obtain the LO resuh it is enough to keep only the 
permutations where the soft gluon e is close by its nearest neighbour, gluon a. 

Note that in the formalism leading to ()25p all the gluons are considered as incoming particles; 
that is, the energies of the gluons c, d, e are negative. In the case when one or two momenta in 
the product {ah) have negative energy, the phase 0^6 is calculated with minus the momenta with 
negative energy, and then n7r/2 is added to (pab where n is the number of negative momenta in 
the spinor product. 

The three jet cross section (jSI) is the square of the matrix element (fT^ calculated using the 
subamphtudes given by (j23). In this way, we obtain 



da = IMP ;^^^'^ n.^, (26) 

where i = a,b, c, d, e and j = c, d, e. To calculate the colhnear LO contribution it is enough 
to keep, in only the permutations where the soft gluon e is the nearest neighbour of 

the incoming gluons a or b. For example, for the case of e colhnear to a we need only retain 
the m(a, e, b, c, d) and m(a, b, c, d, e) subamphtudes, plus the analogous amplitudes with all the 
permutations of the gluons b, c, d. When we sum over the permutations of gluons b, c, d, and 
account for the fact that in colhnear approximation the 4- vector is parallel to a^, we obtain 
the exclusive amplitude of high-ET dijet production. The factor (ae) in the denominator of the 
subamplitude provides the LO logarithm dsae/sae in the cross section. 
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